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ABSTRACT 

Allocation  of  resources  has  become  a  classical 
problem  in  optimization  by  mathematical  programming. 
In  the  field  of  military  applications  attack  aircraft 
assignment  has  been  treated  widely  by  deterministic 
and/or  linear  models.   However,  destruction  of  a  target 
is  no  certainty  nor  is  damage  inflicted  on  a  target 
linear  with  respect  to  the  number  of  weapons  delivered 
on  it.  Recent  extensions  In  the  field  of  nonlinear 
programming  in  conjunction  with  the  widespread  use  of 
electronic  digital  computers  permit  a  more  realistic 
approach  to  this  problem.   This  paper  formulates  a 
stochastic  nonlinear  model  for  assigning  a  force  of 
attack  aircraft  on  a  single  sortie  against  fixed  location 
targets.  The  number  of  aircraft  alive  at  weapon  release 
on  any  pass  of  a  series  against  a  given  target  is  treated 
as  a  random  variable.   The  total  value  of  damage  to  all 
targets  is  taken  as  the  measure  of  effectiveness  and  a 
particular  form  of  the  objective  function  derived.  The 
parameters  of  the  model  and  the  form  of  the  constraint 
equations  are  also  discussed. 
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CHAPTER  I 

THE  PROBLEM 

I.  BACKGROUND 

Allocation  of  resources  has  become  a  classical 
problem  in  optimization  by  mathematical  programming. 
Linear  programming  has  been  used  extensively  to  provide 
a  solution.   It  is  valid  when  the  relationship  between 
variables  is  linear,  however,  linearity  is  not  always 
the  case.   For  example,  an  individual  may  find  the  use- 
fulness derived  from  owning  five  automobiles  not  nec- 
essarily five  times  that  derived  from  owning  Just  one. 
In  the  past,  for  cases  such  as  this,  it  was  customary 
for  one  to  either  be  satisfied  with  a  linear  approxi- 
mation or  to  search  by  other  means  for  a  solution.   The 
extension  of  nonlinear  programming  over  the  past  five 
years  to  its  present  state  provides  the  analyst  with  a 
mathematical  programming  technique  for  treating  the 
allocation  problem  with  certain  types  of  honlinearity 
in  a  direct  rather  than  approximate  fashion.   Dorn  |18| , 
Graves  and  Wolfe  [l?]  ,   Hadley  [25]  „  and  Wolfe  [13]  have 
published  comprehensive  works  on  the  status  of  nonlinear 
programming. 


A  particular  allocation  problem  In  the  field  of 
military  applications  that  has  received  a  great  deal  of 
attention  Is  that  of  attack  aircraft  assignment.   Numerous 
studies  and  papers  have  been  published  on  the  subject, 
and  several  which  have  been  found  to  be  of  particular 
Interest  are!   Limited  War  Operations (U)  |lOJ ,  Tactical 
Air  Warfare  Force  AllocatjLon(U)  [21J  ,  Dresher  U-J  f  and 
Haerlng  [29]  . 

II.  ASSUMPTIONS 

It  was  through  participation  in  the  Center  for  Naval 
Analyses  Study!  Tactical  Air  Warfare,  196*4-.  that  the 
necessity  for  a  formulation  such  as  developed  in  this 
paper  became  apparent.   The  allocation  of  attack  aircraft 
to  fixed  location  targets  is  treated  under  the  following 
broad  assumptions! 

(Al.l)   An  optimal  allocation  is  one  in  which  the 
total  value  of  damage  to  all  targets  is  a  maximum,  subject 
to  the  constraints  imposed  by  the  availability  of  aircraft, 
fuel,  and  weapons. 

(A1.2)   The  damage  inflicted  on  a  target  is  nonlinear 
with  respect  to  the  number  of  aircraft  passes  made 
against  It. 

(A  — ._)  identifies  assumptions^, 
(D  ~.~)  identifies  definitions. 


(A1.3)  There  is  a  probabilistic  interpretation 
of  the  number  of  aircraft  alive  on  a  particular  pass 
against  a  given  target.   That  is,  given  ten  aircraft 
assigned  to  attack  a  target,  there  is  a  probability 
that  all  ten  aircraft  will  not  be  alive  on  any  parti- 
cular pass  against  that  target. 

III.  THE  FORMULATION 

Under  these  assumptions  the  allocation  of  attack 
aircraft  was  formulated  as  a  stochastic  nonlinear 
mathematical  program  with  its  attendant  objective  function 
and  constraint  equations.   The  objective  function  in  this 
case  is  a  nonlinear  function  of  the  number  of  attack  air- 
craft alive  on  each  pass  against  each  target.   Its 
particular  form  is  derived  and  a  discussion  of  the  con- 
straint equations  given  in  Chapter  II . 

The  number  of  passes  against  a  target  is  treated  as 
a  random  variable  in  the  following  way 2   the  number  of 
aircraft  alive  at  release  on  a  particular  pass  against 
a  target  is  a  random  variable,  therefore  the  number  of 
passes  against  a  target  is  a  random  variable  which  is 
the  sum  of  the  number  of  aircraft  alive  at  release  on 
each  pass.   Their  probability  distributions  are  derived 
in  Chapter  III  and  Appendix  A. 

Chapter  IV  is  devoted  to  a  discussion  of  the  parameters 
of  the  model.   These  include  target  value 9  fuel  required, 


racks  available  for  ordnance,  the  probabilities  of 
detection,  engagement,  kill,  acquisition,  and  pass 
survival. 

Conclusions  and  recommendations  comprise  the  final 
chapter  of  this  paper. 


CHAPTER  II 
THE  OBJECTIVE  FUNCTION  AND  CONSTRAINT  EQUATIONS 
I.   THE  OBJECTIVE  FUNCTION 

The  measure  of  effectiveness  chosen  for  the 
allocation  of  attack  aircraft  was  the  total  value  of 
damage  inflicted  upon  enemy  targets.   The  problem  is  to 
maximize  the  total  expected  value  of  damage  subject  to 
certain  constraints 0 

(A2.1)   Assume  that  each  aircraft  assigned  the  i 
target  delivers  the  same  number  of  preferred  weapons  on 
each  pass  against  that  target*  A  preferred  weapon 
implies  that  for  a  specific  target  there  exists  a  weapon 
which  is  most  effective  in  destroying  that  target. 

(A2,2)   On  any  pass  against  a  target  assume  that  the 
target  is  either  killed  or  not  killed.   This  implies  that 
a  target  can  be  destroyed  on  one  pass. 

(A2.3)   Aircraft  make  passes  until  all  weapons  are 
expended  or  until  the  aircraft  is  killed* 

(D2.1)   The  base-aircraft  index  specifies  an  aircraft 
type  located  at  a  specific  base.  For  example, 

j  =  1  denotes  the  set  of  A^  aircraft  from  base 
number  ones 

j  =  2  denotes  the  set  of  F4  aircraft  from  base 


number  one, 


J  =  J  -  1  denotes  the  set  of  A4  aircraft  from 
base  number  S, 

j  ■  J  denotes  the  set  of  F4  aircraft  from  base 
number  S. 

(D2.2)  N.,  denotes  the  number  of  J   base-aircraft 
type  assigned  the  i   target,  i  =  1,2,...,  I  and 

J  =  .L,^,..»,<J. 

(D2.3)   N.  denotes  the  vector  (N,,,  N12>  ...  ,  Njj)» 

(D2.^)  N,  denotes  the  total  number  of  the  j   base- 
aircraft  type  available. 

(D2.5)  Ra*  denotes  the  number  of  passes  per  aircraft 
planned  by  the  J   base-aircraft  type  assigned  the  1 
target, 

(D2.6)  PT.,  ,  a  constant,  denotes  the  probability 
the  i   target  is  killed  by  exactly  one  pass  given  that 
b  preferred  weapons  are  delivered  on  that  pass. 

(D2.7)  W.(N J  denotes  the  random  variable  which 

represents  the  total  number  of  passes  against  the  i 

target.  W. (N. )  will  be  written  as  W.  when  convenient. 

The  range  of  possible  values  is J 

J 
W±   =  0,1,2, ...E  N^Rj,. 

•J  ~"*J» 


(D2„8)  Let  Y.  be  a  random  variable  such  that 
I.  =J0  If  the  1th  target  is  not  killed 

1  if  the  1th  target  is  killed  at  least  once. 


Then 


=  £    (    1   -  PTlb    )w  P  [W1  =  w]    , 


all  w 


and 


p  [i4  - 1]  - 1  -  p  K  =  o] 

-  1  -  £  J  1  -  PT1h)W     P  fw.    -  w]    . 
all  w  1D  L    i  J 

(D2.9)  Let  v.  denote  the  pre-assigned  value  of  the 

th  X 

1   target,,  a*id  let  V  =  £  v.  be  the  total  value  of  the 

i=l  x 
target  complex,, 

(D2.10)  Let  Di  denote  the  random  variable  which 

th 
represents  the  value  of  damage  to  the  i    target  such 

that: 


D^< 

0      If  Y1   s  0 

%    ltjvml 

Then 


M  %A  *•*[»!-*] 


E 

ail  d 


1 

-   v.    (1 


The  total  value  of  expected  damage  (TD)  to  the  target 


complex  Is  I 


(D2.ll)      TD  = 


A E  N 


11  w         r  i 

E     v.    -     E  Y,  E    (   1   -  PT  .     )w     P   [w     =  w] 


i=l      *        1=1   *   all  W 


=  V  -     E  v.  £(1-  PT.,     )W  P  [w.    =     w]    , 

1=1  x     all  w  1D  L    i  J 


w  =   0,1,..,,      E  N. ,R. ,. 
J=l   1J   1J 

To  maximize  TD  one   must  minimize   the   expression 

E     v^     E   (   1  -  PTlb   )w     P  [Wt  =  w]  ,   which 

fch 

is  simply  the  sum  of  the  product  of  the  value  of  the  i 
target  and  the  probability  the  i   target  is  not  killed. 
Therefore  the  objective  function  is 

(2.1)  min  E  v.    E  (  1  -  PT.,  )w  P  [w,  =  wl  . 
1=1  x  all  w        1D      L  i     J 

II.   THE  CONSTRAINT  EQUATIONS 

The  minimization  of  the  objective  function  is  subject 
to  certain  constraints  on  the  resources  available.   The 
four  constraints  given  below  are  a  non-exhaustive  subset 
of  those  which  could  be  imposed,  but  are  representative 
of  the  most  important  ones 3 

1.   The  number  of  aircraft  allocated  to  targets 
must  be  less  than  or  equal  to  the  total  number  of  aircraft 
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available,  I.e., 

(2.2)  jJLl^-P,   .  *U      *  0, 

where  N,  =  the  total  number  of  j   base-aircraft  type 
available,  and  N.  ,  =  the  number  of  J   base-aircraft 

t"  Vi 

type  assigned  the  i   target. 

2.   The  total  fuel  required  must  be  less  than 
or  equal  to  the  total  fuel  available,  i.e., 


(2-3)  &  fuNij '  fj   • 

feVi 

where  F,  ■  the  total  fuel  available  for  the  j   base- 
aircraft  type,  and  f . .  =  the  fuel  required  for  the  j 
base-aircraft  type  to  strike  the  i   target. 

3.   The  total  ordnance  loading  must  be  less  than 
or  equal  to  the  ordnance  available,  i.e., 

(2'*>  &  bUnNlJ  "BJn  > 

where  B.  =  the  number  of  weapons  of  type  n  available  for 
the  J   element,  and  b, .  =  the  number  of  preferred 
weapons  of  type  n  the  j   element  carries  to  the  i 
target.  % 

^.   There  exists  an  upper  limit  on  the  number 
of  aircraft  which  one  is  willing  to  lose  on  a  given 
mission,  i.e., 

(2-5)  Jx  LuNij  -  lj   • 


where  L*  =  the  maximum  acceptable  number  of  J   base- 
aircraft  losses j  and  L.  .   =  the  expected  percentage 
attrition  of  the  J   base-aircraft  type  assigned  the 
1th  target . 

There  are  numerous  ramifications  to  these  constraints 
which  will  not  be  covered  in  detail  here*   For  example 
there  are  additional  restrictions  to  (2.3)9  (2„4-),  and 
to  (2.5)  in  that 9  respectively ,  f. .  is  constrained  by 

the  amount  of  fuel  an,  aircraft  can  carry,  b«  .   is  con- 

ijn 

strained  by  the  maximum  ordnance  load  and  rack  restrict- 
ions of  the  aircraft  type,  and  L,  ,  is  a  function  of  the 
aircraft  type,  speed,  penetration  altitude 9  the  enemy 
defenses,  and  the  number  of  passes  made  against  the 
target. 

The  allocation  problem  has  now  been  formulated  as* 

I 

(2.1)  min  E  v.     £  (  1  -  PT„,  )W  P  fw,  =  w 

1=1  1  all  w        1D      L  * 

subject  to  I 

I 

(2.2)  E  NtJ  =  Nj  9      N^  -  0  , 

I 

(2.3)  %     fijNij  ~  Fj  » 

(2o5)  JxLiAj  *h  . 
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III.   NONLINEAR  PROGRAMS 

Rosen  [?]  ,    111   and  Fiacco  and  McCormick  [20 J   , 
|22]  ,  I  30 1  have  developed  algorithms  which  have  been 
successful  in  solving  nonlinear  programming  problems 
subject  to  linear  or  nonlinear  constraints. 

The  Sequential  Unconstrained  Minimization  Technique 
(SUMT)  of  Fiacco  and  McCormick  and  the  Gradient 
Projection  (GP)  technique  of  Rosen  are  both  available 
through  the  IBM  Share  General  Program  Library  as  Share 
Distribution  3189,  RAC  SUMT,  and  Share  Distribution  1399, 
SD  GP  90  respectively. 

The  use  of  SUMT  is  precluded  in  this  allocation 
problem  due  to  the  non-differentiable  nature  of  the 
objective  function.   That  is,  W.  is  a  discrete  integer- 
valued  random  variable.   The  concluding  remarks  of  the 
Share  write  up  of  RAC  SUMT  Indicate  that  special  sub- 
routines are  being  developed  to  handle  non-differentiable 
functions. 

The  program  GP  90  should  handle  the  non-differentiable 
objective  function  since  it  uses  two-sided  differences 
in  place  of  the  gradients.   Thus  it  is  unnecessary  to 
explicitly  evaluate  the  gradients  of  the  objective 
function. 
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CHAPTER  III 
THE  PROBABILITY  DISTRIBUTIONS 
I.  ASSUMPTIONS 

The  probability  distributions  of  two  of  the  random 
variables,  the  number  of  aircraft  alive  at  release  on 
a  particular  pass  and  the  total  number  of  passes  against 
a  target,  are  presented  in  this  chapter. 

The  following  assumptions  were  made* 

(A3.1)  An  individual  aircraft  is  assigned  only  one 
target  per  sortie* 

(A3. 2)  A  raid  is  composed  of  one  base-aircraft  type, 
but  more  than  one  raid  can  be  assigned  to  a  target, 

(A3. 3)  Enemy  fighters,  if  scrambled  against  a  raid, 
are  sent  in  numbers  sufficient  to  engage  each  aircraft 
in  that  raid. 

II,  DEFINITIONS 

The  following  definitions  were  useful, 

(D3.1)  PD.^n)  denotes  the  probability  that  a  raid 

th 
of  size  n  is  detected  enroute  from  the  j   base-aircraft 

location  to  the  i   target, 

(D3.2)  FE± ,(N,n)  denotes  the  probability  that  a  raid 

th 
of  size  n  is  engaged  enroute  from  the  j  '   base-aircraft 

location  to  the  1   target  given  that  the  raid  is  detected 

12 


and  given  that  a  total  of  N  aircraft  are  employed  in 
strike  operations, 

(D3.3)  PK  denotes  the  probability  that  any  aircraft 
in  a  raid  is  killed  enroute  given  that  the  raid  is 
engaged. 

(D3.*0  PA..(n)  denotes  the  probability  that  a  raid 
from  the  J   base-aircraft  location  finds  the  i   target 
given  that  n  aircraft  survive  enroute. 

(D3«5)  ?Riiv(n)  denotes  the  probability  that  any 
aircraft  in  a  raid  from  the  j   base-aircraft  location 
survives  until  the  k   release  against  the  i   target 

given  that  n  aircraft  are  alive  commencing  the  first  pass, 

st 

k  =  1,  or  given  that  n  aircraft  are  alive  at  the  (K  -  1) 

release  for  the  second  and  subsequent  passes,  K  =  2\ 
k  =  1 ,  c  ,  J  ,  .  •  •  ,xi.  .  • 

(D3»6)  Xi1k(N. m)   denotes  the  random  variable  which 
represents  the  number  of  the  J   base-aircraft  type  alive 
at  release  on  the  k   pass  against  the  i   target  given 
that  N  .  are  assigned.   The  range  of  possible  values  is! 
X1Jk(N1;J)  =  0,1,2,...,N1;J.   For  simplicity  X1Jk(N1;J)  is 
sometimes  written  as  X..,. 

III.  DISTRIBUTION  OF  THE  NUMBER  OF  AIRCRAFT 
ALIVE  AT  RELEASE  ON  THE  FIRST  PASS 

Under  the  above  assumptions  the  event  that  x  aircraft, 
i»0,  are  alive  at  release  on  the  first  pass  against  a 
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target  can  occur  In  three  mutually  exclusive  and 
exhaustive  ways" 

A  raid  of  size  N  is  undetected  enroute,  finds  the 
target  and  x  aircraft  among  N  survive  target  defenses 
on  the  first  pass, 

OR  a  raid  of  size  N  is  detected  enroute9  but 
unengaged,  finds  the  target  and  x  aircraft  among  N 
survive  target  defenses  on  the  first  pass, 

OR  a  raid  of  size  N  is  detected  enroute9  engaged, 
some  n  aircraft  survive  enroute,  n  -  x,  those  n  find  the 
target  and  then  x  among  n  survive  target  defenses  on  the 
first  pass. 

The  extension  of  this  type  of  reasoning  used  in 
conjunction  with  the  laws  of  elementary  probability 
theory  produced  the  probability  distribution  of  the 
number  of  aircraft  alive  at  release  on  the  first  pass, 
X-.,,  presented  in  Table  I. 

IV.  DISTRIBUTION  OF  THE  NUMBER  OF  AIRCRAFT 
ALIVE  AT  RELEASE  ON  SUBSEQUENT  PASSES 

The  probability  distribution  of  X.1k9  k  ■  2 ,  was 
obtained  from  the  distribution  of  X. ._  by  making  the 
following  observations^ 

The  stochastic  process,  |  X^^  »,  k  ■  19293,...,R .  .7 

is  a  discrete  parameter  finite  Markov  chain.   That  is, 

P  LXiJk  ::  xk  I  Xtjl  ::  xl8  Xij2  =  x2  5  00°  9  Xijk-1  =  xk-l] 

=  P  LXijk  =  xk  I  Xijk-1  =  xk-lJ 

Ik 


The  number  of  the  J   base-aircraft  type  alive  at 

release  on  the  k   pass  against  the  1   target  Is  less 

st 
than  or  at  most  equal  to  the  number  alive  on  the  (k-1) 


kijk  "  Xijk-1      k  "  2»3,^,...,R1J 


P    Kjk  "  *]      '  j^      P   [Xijk  =  *  |  Xijk-1  "  m]   P  [Xijk-1  "  m] 


pass.   That  is, 

Therefore  the  conditional  distribution  of  X-..,  k  -  2, 
given  that  X.  ..  -,=m,  is  binomial  with  parameters 
PR..k(m)  and  m  as  presented  in  Table  II. 
By  the  theorem  of  total  probabilities  I 

r 

m=x 

The  synthesis  of  these  facts  yielded  the  probability 
distribution  of  Xi1k,  k  3*  2,  presented  in  Table  III. 

V.   DISTRIBUTION  OF  THE  TOTAL  NUMBER 
OF  PASSES  AGAINST  A  TARGET 

The  derivation  of  the  probability  distribution  of 
W.  necessitated  one  further  assumption! 

(A3«5)  The  number  of  a  given  base-aircraft  type  alive 
at  release  on  any  pass  over  a  given  target  is  statistically 
independent  of  any  other  base-aircraft  types  alive  over 
the  same  target  on  any  pass.   That  is,  X.  ..  is  assumed 
independent  of  X. ./,/ provided  J  f   j; . 

Under  all  foregoing  assumptions  the  probability 
distribution  of  W.  as  presented  in  Table  IV  is  derived 
in  Appendix  A. 
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CHAPTER  IV 
THE  PARAMETERS 

I.  TARGET  VALUE 

The  value  of  a  target  may  be  characterized  by  a  time 
dependent  indexg  ioe.?  target  value  will  change  as  the 
tactical  situation  changes  over  time.   Furthermore ,  the 
relative  values  of  targets  in  given  tactical  situations 
will  depend  strongly  upon  the  doctrine  adopted  by  the 
local  military  commander0   If  an  index  of  target  value 
is  to  provide  a  basis  for  decisions,  it  must  be  formulated 
with  consideration  given  to  the  whole  social;,  political, 
economic^  and  ethical  processes  within  which  military 
action  takes  place c   The  intent  of  this  section  is  not 
to  formulate  a  scheme  for  target  value  indexing  but  to 
present  some  aspects  of  the  problem  which  should  be 
considered  in  any  indexing  plan0 

Kaplan  |2|  delimits  areas  of  empirical  study  by 
which  each  element  in  the  concept  of  military  worth  can 
be  measured  or  appraised,.   These  areas  include  values , 
objectives ,  welfare 9  achievement 9   and  wortho 

Kaysen  |3|  outlines  criteria  for  selecting  target 
systems  as„ 


20 


1,  Military  importance  --  a  rough  classification 
of  the  value  to  enemy  military  operations  of  all  types 
of  equipment  and  supplies  used  by  enemy  forces  J 

2#  Percent  direct  (or  indirect)  military  use  —  the 
share  of  total  output  of  a  product  or  service  which  goes 
into  military  use  J 

3.   Depth  —  an  indication  of  the  time  available  to 
the  enemy  for  the  organization  of  substitute  consumption, 
alternate  production,  etc.,  before  suffering  military 
damage \ 

k*     Economic  vulnerability  —  includes 

a.  Ratio  of  capacity  to  output, 

b.  Substitutability  for  processes  and  equipment, 

c.  Substitutability  for  product, 

d.  Process  and  plant  layout  vulnerability, 

e.  RecuperabilityJ 

5.   Physical  vulnerability; 

6#   Location  and  size  of  the  target  system. 

In  addition  to  the  above  criteria  for  indexing  target 
value  the  effect  on  enemy  morale  should  be  considered. 

Hesse  and  Mitchell  in  "Limited  War  Campaigns! 
Dacca  Method' •  |  10 J  have  derived  an  index  of  target  value 
which  incorporates  several  of  the  criteria  mentioned 
above.   They  consider  military  targets  from  a  balanced 
ton  standpoint  using  weighting  factors  for  consumable  and 
reuseable  equipment,  supplies,  and  manpower.   This  index 
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is  quite  involved;,  but  it  provides  a  method  of  assigning 
target  values  for  most  conceivable  military  targets „   It 
should  serve  well  as  a  useful  working  tool  for  a  first 
look  at  target  value  in  the  allocation  problem,, 

II o  FUEL  AND  ORDNANCE  LOADINGS 

To  maximize  damage  to  the  enemy  requires  that  air- 
craft be  put  over  their  targets  with  optimal  ordnance 
loadings.   A  basic  prerequisite  to  attainment  of  this 
objective  is  determination  of  fuel  required  to  carry  out 
the  given  mission,, 

The  following  basic  mission  profile  has  been  postulated 
to  determine  fuel  requirements, 

1.  Warm-up  and  take-off 8 

2.  Climb  to  cruise  altitude  and  cruise  to  descent 
point  for  run-in  to  the  target 8 

3.  Descent  to  target  run-in  altitude 9 
bo      Run-in  to  the  target 9 

5o  Ordnance  delivery , 

6.  Run-out  from  the  target D 

7.  Climb  to  cruise  altitude  and  return  to  base9 

8.  Descent  to  the  base  for  landings 

9.  Landing  and  reserve „ 

For  ease  of  computation  the  following  assumptions 
have  been  made" 

(A40l)   Aircraft  carry  full  internal  fuel  and  use 


fully  loaded  300  gallon  external  fuel  tanks  as  required* 

(A^.2)   A  fuel  weight  of  six  and  one-half  pounds 
per  gallon  was  used; 

(A^.3)   No  fuel  is  burned  and  no  distance  over  the 
ground  is  covered  in  descents; 

(A*J-.^)   All  external  store  racks  for  ordnance  are 
equivalent  on  any  particular  aircraft  type  J 

(A4--5)   There  are  no  asymmetric  load  restrictions 
for  launching  aircraft; 

(A4-6)   Five  minutes  fuel  at  normal  rated  power  was 
allotted  for  warm-up  and  take-off; 

(A4--7)   Thirty  minutes  fuel  at  sea  level  maximum 
endurance  was  allotted  for  landing  and  reserve. 

Utilizing  the  mission  profile  and  the  assumptions 
stated  above  one  may  formulate  the  basic  fuel  required 
equation  as 

(**.l)   Total  fuel  required  =  Take-off  fuel  +  (2  x 

cruise  fuel  rate  x 
(Range  -  Climb  distance  - 
Run-in  distance))  +  (2  x 
Run-in  distance  x  Run-in 
fuel  rate)  +  (Time  over 
target  x  Target  fuel  rate) 
+  Landing  and  reserve  fuel, 
After  determining  the  total  fuel  required  the  next 
step  is  to  ascertain  if  the  fuel  required  is  less  than 
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the  fuel  available.   Successful  mission  performance  is 
constrained  by 

(4.2)  Fuel  available  -  Fuel  required  ^0. 
Substituting  for  Fuel  available  in  (4.2)  one  obtains 

(4.3)  (Internal  fuel  +  External  fuel)  -  Fuel 
required  ^0. 

Solving  for  External  fuel  results  in 

(4.4)  External  fuel  »  Fuel  required  -  Internal  fuel. 
Once  (4.4)  has  been  solved  it  is  possible  to  determine 

the  number  and  size  of  external  fuel  tanks  required  and 
subsequently  to  determine  the  number  of  external  store 
stations  available  for  ordnance  loading. 

The  maximum  ordnance  load  may  be  obtained  by  solving 

(4.5)  Ordnance  load  =  Maximum  take-off  weight  - 

(Basic  aircraft  weight  + 
Internal  fuel  +  External  fuel). 
Appendix  B  contains  a  program  (PROGRAM  0RDL0AD)  written 
in  FORTRAN  IV  for  the  CDC-1604  digital  computer  which 
computes  for  a  given  aircraft  mission  the  fuel  required, 
the  number  of  store  stations  available  for  ordnance,  and 
the  maximum  ordnance  loading.   This  program  uses  the 
assumptions  and  the  algorithms  outlined  above  for 
computations.   Inputs  required  to  the  program,  definitions, 
a  program  listing,  and  a  sample  output  are  given  in 
Appendix  B. 
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Ill,  OTHER  PARAMETERS 

The  probability  that  an  aircraft  survives  to  make 
passes  against  a  target  is  a  function  of  the  probability 
of  detection  by  the  enemy 9  the  probability  of  engagement 
by  the  enemy,,  the  probability  of  being  killed  enroute  to 
the  target  by  the  engaging  force „  the  probability  of    -  . 

acquiring  the  target 9  and  the  probability  of  being  alive 

th 
at  the  weapons  release  point  on  the  k   pass  against  the 

targe to   Values  for  all  of  the  above  are  required  to 

obtain  the  distribution  of  W,  0   To  minimize  the  objective 

function  one  also  needs  values  for  the  probability  that 

the  i   target  is  killed  by  one  pass  given  tVnt  b 

weapons  are  released  on  that  passc   These  probabilities 

will  be  considered  as  input  parameters 0 

Derivation  of  each  of  these  probability  distributions 
is  beyond  the  scope  of  this  papery  however 9  the  remainder 
of  this  section  will  be  devoted  to  discussion  and/or 
references  to  derivations  of  these  distributions,, 
Throughout  the  investigation  of  these  parameters 9  past 
data  will  in  some  instances  be  adequate  and  reliable,  and 
the  parameters  could  be  estimated  by  Bayesian  techniques. 
World  War  II „  Korea 9  and  Viet  Nam  experiences  should 
result  in  extensive  information  which  may  be  used  to 
increase  confidence  in  the  estimates  used. 

Koopman  1  derives  an  expression  for  detection 
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probability  both  for  visual  and  for  radar  detection. 
Electronic  and  other  intelligence  data  coupled  with 
characteristics  of  friendly  r^dar  should  allow  an 
estimate  to  be  obtained  for  detection  probability.   Crout, 
Fay,  and  Harvey  pi  have  formulated  a  model  describing 
a  bomber's  penetration  into  hostile  territory  in  the 
face  of  a  given  area  defense  strategy.   Specifically  an 
expression,  is  developed  to  compute  the  probability  that 
an  attacking  aircraft  is  killed  before  penetrating  to  a 
given  depth.   This  model  then  accounts  for  the  probability 
of  detection,  engagement,  and  kill. 

Simultaneously  with  passing  through  enemy  area 
defenses  friendly  aircraft  must  attempt  to  acquire  their 
assigned  targets.   Acquisition  includes  the  processes 
of  search,  detection,  identification,  and  flying  the 
aircraft  into  a  position  to  make  a  weapons  pass  against 
the  target.   Visual  acquisition  of  the  target  is,  of 
course,  the  most  reliable  and  informative  method  of 
locating  a  given  target.   Due  to  circumstances  such  as 
weather,  terrain,  camouflage,  aircraft  speed,  and  aircraft 
altitude  a  pilot  will  sometimes  have  to  rely  on  radar, 
infrared,  microwave,  or  vectoring  from  personnel  not  in 
the  aircraft  to  acquire  his  target.  A  combination  of 
these  methods  should  increase  the  probability  an  aircraft 
acquires  its  target.   Erickson  [9]  presents  a  study  on 
the  visual  capabilities  of  a  pilot  searching  for  ground 
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targets,   Kamrass  and  Heckroth  |15J  include  a  section 
on  target  acquisition  which  discusses  detection 
probabilities  and  an  analysis  of  types  of  sensors  and 
methods  of  operation  which  will  be  most  effective  for 
reconnaissance  and  strike:  missions^ 

Once  an  aircraft  has  acquired  a  target  and  is  positioned 
for  a  weapons  release  pass  it  will  usually  be  in  the 
area  of  effectiveness  of  the  point  defenses  of  the  target. 
These  point  defenses  can  run  the  gamut  from  no  defense 
to  a  completely  automatic  system  of  surface-to-air 
missiles  with  a  ground  level  to  outer  space  capability,, 
The  probability  an  aircraft  can  successfully  release 
weapons  over  the  target  is  highly  dependent  upon  the  type 
of  defensive  weapon  encountered „   Selection  of  Aircraft 
for  Tactical  Air  Missions   2^  discusses  and  gives  values 
for  probabilities  of  kill  and  the  effect  of  various 
penetration  aids  associated  with  several  types  of  point 
defense  weapons 0 

If  an  aircraft  survives  the  point  defenses 9  it  will 
release  weapons  against  the  target <>   The  probability  of 
killing  the  target  on  one  pass  will  be  a  function  of  the 
type  and  number  of  weapons  released  on  that  pass„  the 
skill  of  the  individual  pilot,,  the  weather 9  the  type  of 
terrain9  and  the  type  of  target c   The  Naval  Ordnance  Test 
Station  (NOTS)  has  worked  on  the  development  of  weapons 
and  the  formulation  of  models  which  allow  the  probabilities 


of  kill  of  those  weapons  to  be  computed.   NOTS  technical 
publications  include  those  of  a  probabilistic  nature  — 
Kusterer  [a]  ,  W.B.  Simecka  [l6J   ,  Verry  [19]  ,  K.D. 
Simecka  [26   ,  and  Strang  I  28  J  J  and  those  of  a 
computational  nature  —  Weldon  and  Young  |  14   ,  W.  B. 
Simecka,  et.  al.  [23]  »  and  ''Conventional  Air-Delivered 
Strike  Weapons"  [27]  . 

An  alternative  to  derivation  of  the  various 
probability  distributions  for  all  the  parameters  would 
be  to  use  three  estimates  of  each  parameter,  i.e.,  a  best, 
an  optimistic,  and  a  pessimistic  estimate.   These  will 
be  referred  to  as  BOP  estimates.  After  obtaining 
solutions  to  the  allocation  problem  using  these  BOP 
estimates,  a  sensitivity  analysis  should  be  performed 
to  discover  which  parameters  significantly  affect  the 
results.   Following  this  analysis  high  confidence 
estimates  for  the  most  sensitive  ''parameters''  could 
be  obtained  using  the  distribution  theory. 


28 


CHAPTER  V 
CONCLUSIONS  AND  RECOMMENDATIONS 
I.  CONCLUSIONS 

This  analysis  has  demonstrated  that  the  allocation 
of  aircraft  to  fixed  location  targets  is  a  complex 
problem,.   If  it  had  been  reasonable  to  assume  that 
attack  aircraft  make  a  single  pass9  as  in  strategic 
nuclear  warfare 9  then  the  nonlinear  model  reported  by 
Fiacco  and  McCormick  |22J  could  have  been  used,,   However 
in  limited  warfare 9  for  which  this  model  was  developed,, 
it  is  not  reasonable  to  assume  a  single  passG   Furthermore, 
we  believe  that  a  deterministic  treatment  of  repeated 
passes  would  not  suitably  reflect  the  expected  outcome 
in  even  an  approximate  fashion,.   The  substantiation  of 
this  assertion  is  dependent  upon  the  generation  of  the 
probability  distribution  of  W., 

The  objective  function  for  this  model  is  non- 
dif ferentiable"  this  property  restricts  the  class  of 
usable  and  currently  available  nonlinear  programming 
algorithms  to  that  of  Rosen  [?]  „ 
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II.  RECOMMENDATIONS 

It  is  recommended  that . 

1.  This  model  be  coded  for  use  on  an  electronic 
digital  computer  utilizing  the  algorithm  of  Rosen  I7I  • 

2.  The  list  of  targets  be  grouped  Into  type  or 
character  classes  and  an  Investigation  made  of  the  feas- 
ibility of  using  common  values  within  each  group  for  the 
parameters  PA,.(n),  PR.  .,  (n)  and  v.  respectively. 

3.  The  list  of  targets  be  grouped  Into  geographical 
classes  and  an  Investigation  made  of  the  feasibility  of 
using  common  values  within  each  group  for  the  parameters 
PD.  An)   and  PE..(N,n)   respectively. 

^.   A  sensitivity  analysis  using  BOP  (best, 
optimistic,  pessimistic)  estimates  of  the  parameters  be 
made  to  de  term  Jme  those  requiring  further  analysis. 

5.   The  feasibility  Of  invoking  the  Central  Limit 

Theorem  on  the  expression  W,  =  £  E  X.  -.  be  determined. 

1    J  k  1JK 
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APPENDIXES 


APPENDIX  A 
DERIVATION  OF  THE  PROBABILITY 
DISTRIBUTION  OP  W^N^) 


,  or 


The  event  W.  =  0  occurs  if  and  only  if  £  %   X- *k  ■  0 
if  and  only  if  X.  ..  =  0  for  all  J  and  all  k.   It  was  noted 
that  X1Jk  ^  x1Jk-1»  k  ~   2,3,^,...,R1;J,  so  that  Xljk  =  0 
for  all  j  and  all  k  if  and  only  if  X^  =  0  for  all  J. 
It  was  assumed  that  X.  .,  is  Independent  of  X.  ./,  /  provided 
that  J  4   y  >  so  that: 

P  [Wi  =  o]  =  p[s  £  ^jk^^^ijk3  °  for  a11  J  and  a11  k] 
J  k 

-  P[xln  =  0  for  all  i]   =  ,  p[xln  =  o] 

The  distribution  of  W.  ^  1  was  obtained  in  a 

slightly  different  manner.   The  event  W.  =  1  occurs  if 

and  only  if  for  exactly  one  J,  X.  .,  =  1  and  X.  12  -  0» 

but  the  event  W.=  2  occurs  if  for  exactly  one  J, 

X-.,  =  2  and  X-  12  =  °»  or  if  for  exactly  one  j,  X..,  =  1»X.  .g8^ 

and  X.  .-  =  0,  or  if  for  exactly  two  J's,  X.  ,,  =  1  and 

X. 12  =  0«   This  enumeration  can  be  continued  and  was  up 

to  W.  =  4,  it  produced  a  lengthy  table  from  which  a 

recursion  relation  was  developed.   The  use  of  this  relation 

in  conjunction  with  the  probability  distribution  of  X.,, 

produced  the  probability  distribution  of  W  in  the  following 

manner. 
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(D6.1)   Let  faco  denote,  for  all  combinations  of.   Where 
• 'combinations  * s  is  used  in  the  sense  of  combinatorial 
analysis,  e.g.,  Parzen  16  9    then 

p  [wt  =  i]  =  e[p  xljX  =  i]  p[xij2  =  o  |xin  -  i]^    p[xlnl  =  o] 


p[wi =  2J =  *  .  K  p[xui =  2-aJ  p[xu2  =  ai  i  xui  = 2  -  aJ 


J   a1=0 


.p[x, 


J3 


=  0 


Xij2   =  aJ      »,   P[Xinl  =  °] 


n^j 


\JL  p  ixi).i  =  x]  p  K^z  =  °  I  xua  "  i] 


faco 


■    -P[XiJal  =1]    P   [XIJ82   -   0  Kj.l  "   X] JJ:       ,      P[Xinl  =    °] 


NOTE 


:      P[Xijk  =   Xl  Xijk-1  =  m]    = 


n^jig Ja 
0  for  all  x  >  m  and 


P[X13k  =  °  I    Wl  =  °]     -1.   Table  "• 


2      a,-l 


P  t"l  "  31   -  ?  a  A  .!*o     '  k»  =  3   "  ^ 


1    v  "2 


.P[X1J2   =ax  .  a2  |  Xin  =  3    -  ajp[xij3   =  a£  |  Xlj2   =  ^   -  aJ 
•P[Xij^  =0  I  XiJ3   =  a2-l      3      Pt 


Snl  =  °] 


(continued  on  the  following  page 
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2        3-an-l     a„-l 


+  E 

faco     a. 


E  E1  XE        P  [X  =   3    -  a      -  b      J 

=1  bn-o      b21=o  U*1 


•  P[XiJx2   =  bll  I  Xija  =3    -  a,   -  ^Jp[xi3i3  =   0 

•  P[XiJ8l   =  al    "   b2lJ    P[XiJa2    =  b21  I  Xljal   =  al    *   b2lJ 

■ 

•  P[X1J.1  =  1]      P[X1J,2=   0.1    W1] 


XiJ»2   *  bllJ 


•      P[XU3l  =  1]    P[X1J32   "   0    I    XlJ3l  =  1L  1 


nj^Ji,  J  a  i  J  3 

* 


P[Xinl  =  0] 
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By  induction  the  general  term  of  the  probability 
distribution  of  W.  was  obtained  and  is  presented  in  Table 
IV,  where  all  sums  are  defined  in  the  positive  sense, 


i.e.,  the  sum  £  f(y)  is  defined  if  and  only  if  e  ^  d 
y=d 


and  the  following  definitions  were  used! 
(D6.2)   R  denotes  the  max  (R«  J 


l.J 


U 


(D6.3)   r  denotes  the  max  (R9J 
Summation  Operators, 


(D6.*0      B 


i 


b  -1 

r~l  Dik„-1  x 


TT 

k1=2 


ik, 


0 


(D6.5)      C^     =     tt 


TT 

k=2 


w=l 
(D6.6)      A,   «       2 
a  1*0 


Lk-i; 


ak=j-k 


r-1 

n 
k=2 


a        -1 

ak-V 

ak=0 


w-1 
(D6.7)      A2   =       E 

al=sl 


w-a, -1 

XE        B, 

bll=0 


a,-l 

XS        B, 

b21=o 


w-1 

(D6.8)      A,   ■        Z 


J 


a^j-1 


j-1  a      -.-a.-l 

n        t=1E   t        B. 

t-i    btl=o 


3   -  3 B^s« . • 9 J  and  aQ  -  w 


aUl"1 
bU=0 


J 
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Summation  Operands ', 

(D6.9)     GtB(i)   =  P  [X,         =  pM  -  at(X-dtm)    -  btl] 


[Xi3t2   " 


btl   "  bt2    |    hitX  =  at-l   -  at(1-dtn.)    "  btl] 


r-2 


kt-2 


:,  pKj 


tV1  "  btkt  "  b*kt+1  !  ^tV  "^t-1  "  btki 


p  [x1. 


Jtr        "tr-1  I    "lJtr-l        "tr-2        "tr 


J 


Lxijtr+1 


=   0        X 


ij.r       btr-lJ 


where   d. ,    =  <    0      if   i 


1     if  i  =   J 


r-2 


(D6.10)      Pn     =  £   2±     P  Kjxk  -  ak-l   "  aJ 
P  L*U»*+1  =  ak  ~  ak+l         Xijxk  =  ak-l   "  aki 
P  LXijAr  =  ar-l       XiJir-l  =  ar-2   "  ar-lJ 


P  [XiJxr+l 
(D6.ll)     P 


=    0 


XiJ>r  «  ar 


-J        "  P  tXinl  =   °] 


n^j 


im 


£ 
faco 


I  G*»(1) 


TT 


n#jx , . .. , 


p  [xmi  "  °] 


m  =  2^3j;..jJ-1  and  J  ^  3.' 


(D6:i2)      T±j  -       it       GtJ(i)  J  -  2 

u=JL 
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Therefore  under  the  above  definitions ,  given  values 
for  R  and  J,  the  probability  distribution  of  W.  presented 
In  Table  IV  can  be  expressed  In  closed  form  for 
computational  purposes 0 
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APPENDIX  B 
PROGRAM  ORDLOAD 

The  FORTRAN  IV  program  to  compute  the  fuel  and 
ordnance  loadings  is  presented  in  this  Appendix.   A 
brief  description  of  the  imput  data,  the  definition  of 
variables,  a  program  listing,  and  a  sample  output  will 
be  given.  Data  for  the  sample  problem  was  selected  at 
random,  and  the  resultant  output  bears  little  semblance 
to  any  realistic  situation  that  might  occur.   The  program 
has  been  dimensioned  to  handle  six  aircraft  type,  six 
bases,  and  fifty  targets.   These  dimensions  may  be 
expanded  as  necessary  to  suit  a  particular  user's  need. 
Common  blocks  have  also  been  established  to  facilitate 
Incorporation  of  this  program  as  a  sub-routine  in  a 
larger  allocation  problem. 

The  data  input  is  presented  below.   Variables  with 
subscript  (I)  require  one  data  cefrd  for  each  aircraft 
type. 

Card  i:   FORMAT  (3H0,F10.  0) : 

Fields  l-io:   NTYPE  =  the  number  of  different 
aircraft  type,  1  *  NTYPE  *  6; 

Fields  11-2 O:   JBASE  =  the  number  of  friendly 
bases,  1  ^  JBASE  ±   6; 
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Fields  21-302   KTGT  f  the  number  of  targets , 
1  ^  KTGT  ^  50; 

Fields  31-402   TIME  p  the  time  (  tenths  of 
hours  )  over  the  target  for  delivery  of  ordnance. 
Card  28.<>.!,NTYPE+12   FORMAT  (  5F10„0SI10  )2 

Fields  1-102   GAXTOW(I)  =  maximum  gross  aircraft 
take-off  weight  (pounds); 

Fields  ll-20:   BACFTW(I)  =  basic  aircraft  weight 
(pounds) I 

Fields  21-302   ENTFW(I)  =  maximum  internal  fuel 
capacity  (pounds) ; 

Fields  31-402   EXTFW(I)  =  maximum  external  fuel 
capacity  (pounds) ; 

Fields  41-502   GAXBMW(J)  =  maximum  gross  ordnance 
load  (pounds); 

Fields  51-602   IRACK(I)  -  maximum  number  of 
external  store  racks  on  an  aircraft; 

Card  NTYPE+25).„o9(2xNTYPE)+12   FORMAT  (5F10„0)2 

Fields  l-io:   FTAKE(I)  =  fuel  (pounds)  required 
for  warm-up  and  take-off; 

Fields  11-202   FCLIMB(I)  ■  fuel  (pounds) 
required  for  climb  to  cruise  altitude; 

Fields  21-302   FLOIT(l)  =  loiter  fuel  rate 
(pounds  per  minute )—  not  required  in  the  current 
formulation  of  the  program; 
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Fields  31-1+0:   DASHSP(I)  =  speed  (knots)  over 
the  target  at  military  power; 

Fields  ^1-50:   FLAND(I)  =  landing  and  reserve 
fuel  (pounds). 

Card  (2xNTYPE)+2,...,(3xNTYPE)+i:   FORMAT  (2F10.0): 

Fields  l-io:   DISTML(I)  =  distance  (nautical 
miles)  covered  in  the  climb  to  cruising  altitude. 

Fields  ll-20:   DISTRI(I)  =  run-in  distance 
(nautical  miles)  to  the  target. 

Card  (3NTYPE)+2:  FORMAT  (  2F10.0): 

Fields  l-io:  ALT1  =  cruise  altitude  (thousands 
of  feet); 

Fields  ll-2  0:  ALT2  =  run-in  altitude  (thousands 
of  feet). 

Card  (3xNTYPE)+3,...,(4xNTYPE)+2:   FORMAT  (^FIO.O): 

Fields  1-10:   SFC(I).,=  sea  level  cruise 
specific  fuel  consumption  (pounds  per  nautical  mile  = 
ppnm); 

Fields  ll-20:   SFCMIL(I)  =  sea  level  military 
power  specific  fuel  consumption  (ppnm); 

Fields  21-30:   SFCK(I)  =  cruise  specific  fuel 

consumption  altitude  correction  factor  (ppnm  per  thousand 

feet  altitude) ; 

Fields  31-^0:  SFCMILK(I)  =  military  power 

specific  fuel  consumption  altitude  correction  factor 

(ppnm  per  thousand  feet  altitude). 

k2 


Card  (^xNTYPE)+3P..o9 (^xNTYPE 3+2+ (JBASExKTGT):: FORMAT  (F10.0): 

Field  1-10S   RANGE (J9K)  =  range  (nautical  miles) 
from  Jth  base  to  the  Kth  target. 

Definitions  of  other  variables  used  In  the  program 
areT 

FCRUSA(I)  =  cruise  altitude  fuel  consumption  rate 
(ppnm); 

FCRUSL(I)  =  run-in  altitude  fuel  consumption  rate 
(ppnm); 

FDASHM(I)  ■  run-in  altitude  military  power  fuel 
consumption  rate  (ppnm) I 

FAVAIL(I)  =  maximum  total  fuel  available  for 
aircraft  type  i; 

FRQRD(I9J9K)  =  fuel  required  for  aircraft  type  I 
to  complete  mission  from  base  J  to  target  K  and  return; 

EXTFRQ(I5J?K)  =  external  fuel  required  for  air- 
craft type  to  complete  mission  from  base  J  to  target  K 
and  return^ 

KRACK(I9 J9K)  =  the  number  of  external  store  racks 
available  on  aircraft  type  I  to  target  K  from  base  j; 

ORDLD(I9J9K)  -  ordnance  load  (pounds)  aircraft 
I  can  carry  to  target  K  from  base  J. 

A  complete  program  listing  and  a  sample  output 9  which 
uses  the  data  at  the  end  of  the  program  listing,  are  given 
on  the    last    four  pages  of  this  Appendix,,   The  sample 
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output  tables  the  aircraft  number,  the  base  number, 
the  target  number,  the  range  from  base  J  to  target  K, 
the  total  fuel  required  (lbs.),  the  Internal  and  external 
fuel  required  (lbs.),  the  number  of  ordnance  racks  available, 
and  the  ordnance  load  capability  (pounds). 
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